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Abstract  -  The  problem  of  the  scattering  of  a  continuous 
acoustic  wave  (CW)  from  a  submerged  object  in  a  waveguide 
has  received  extensive  treatment  in  the  frequency  domain, 
where  phenomena  associated  with  the  resonant  behavior  of 
the  object  have  been  studied  in  detail.  This  paper  concentrates 
on  the  time-domain  interpretation  of  the  scattering  problem. 
Computation  of  the  scattered  field  from  an  incident  pulse, 
which  requires  calculations  for  many  frequencies,  is  performed 
using  a  new  perturbation  technique  (computed  to  all  orders) 
that  allows  rapid  calculation  of  the  waveguide  modes.  An 
implementation  of  the  inverse  Fourier  transform  is  presented 
that  is  more  appropriate  to  resonance-type  spectra.  The 
computational  tools  developed  are  used  to  elucidate  various 
physical  aspects  of  pulse  scattering  in  a  waveguide,  including: 
group  velocity  behavior,  direct  modal  contributions,  and 
resonance  contributions  from  individually  excited  modes. 

1.  INTRODUCTION 

The  frequency  domain  treatment  of  acoustic  scattering  from 
objects  in  a  bounded  medium  such  as  an  ocean  waveguide 
has  received  extensive  study  [1].  It  is  our  purpose  in  this 
paper  to  present  results  on  extending  this  work  to  the  time 
domain.  The  essential  steps  in  the  process  arc  to  1)  calculate 
the  single-frequency  incident  pressure  field,  2)  compute 
the  scattered  field  off  an  object  of  interest  and  3)  repeat  the 
calculations  for  many  frequencies  and  synthesize  the  scattered 
field  by  Fourier  integrating  the  results  into  a  time  domain 
response.  This  is  physically  equivalent  to  transmitting  and 
scattering  a  pulse  off  the  target. 

To  accomplish  the  above,  we  compute  the  incident  pressure 
field  in  the  ocean  waveguide  using  a  new  perturbative  approach 
to  normal  modes  [2].  The  modes  provide  a  set  of  orthogonal 
basis  functions  in  which  to  expand  the  incident 
single-frequency  field.  We  then  use  the  T-Matrix  method  to 
couple  the  incident  field  to  the  scattered  field.  If  this  calculation 
is  repeated  for  a  number  of  frequencies,  an  inverse  Fourier 
transformation  of  the  complex  pressure  field  will  produce  a 
time  domain  scattered  field.  By  reconstructing  only  certain 
combinations  of  modes,  or  individual  modes,  we  can  break 
up  the  time  domain  response  into  various  modal  contributions. 
We  can  also  produce  pulses  for  various  frequency  ranges. 
The  ability  to  manipulate  the  frequencies  and  modes  allow  us 
to  decompose  the  time  domain  response  (which  can  be  very 
complicated)  into  its  constituents,  thereby  elucidating  the 


physical  mechanisms  of  the  different  arrivals.  This  should  be 
especially  useful  in  the  case  of  resonance  scattering  from 
elastic  targets. 

Previous  work  has  been  done  on  the  time  domain  calculation 
of  scattering  from  unbounded  elastic  objects  [3].  Our  goal  is 
to  extend  this  to  elastic  objects  in  a  waveguide.  Since  we 
already  have  the  T-Matrix  codes  necessary  to  compute  the 
scattered  expansion  coefficients,  the  additional  tools  we  re¬ 
quire  are  a  way  to  calculate  normal  modes  rapidly,  yet 
accurately,  and  a  numerical  method  for  accurately  evaluating 
the  inverse  Fourier  transform  when  resonance  effects  are 
important.  In  this  paper  we  will  report  on  progress  we  have 
made  in  developing  these  tools,  with  primary  emphasis  on 
the  normal  mode  calculations.  We  have  developed  a  new 
perturbation  method  that  allows  rapid,  accurate  calculation  of 
the  normal  modes  from  which  the  complex  pressure  field  can 
be  computed  [2].  The  new  perturbation  method  allows  us  to 
write  the  perturbed  normal  modes  as  a  linear  combination  of 
sine  functions,  which  is  important  in  allowing  us  to  expand 
the  incident  field  into  a  spherical  representation  for  coupling 
to  the  T-matrix. 

In  the  following  sections  we  will  present  the  basics  of  the 
T-Matrix  approach  to  scattering,  the  extension  of  this  technique 
to  the  time  domain,  the  basic  ideas  behind  the  new  perturbative 
approach  to  normal  modes,  a  new  frequency  integration 
method,  and  some  first-order  results  illusuating  the  kinds  of 
quantities  we  can  compute. 

2.  EXTENDING  THE  T-MATRIX  TO  TIME  DOMAIN 

The  basic  reference  for  the  following  material  is  Werby 

[4] .  In  this  section  we  will  review  the  methods  by  which  the 
waveguide  pressure  field  is  coupled  to  the  scattered  near  field 
from  the  target.  The  Extended  Boundary  Condition  method 

[5]  (or  T-mauix)  transforms  the  incident  field  into  the 
scattered  near  field.  The  incident  field  is  expanded  into  a 
spherical  representation,  where  the  expansion  coefficients  are 
represented  by  the  column  vector  a.  The  T-matrix  transforms  a 
into  f,  the  expansion  coefficients  for  the  scattered  near  field,  or; 

f  =  Ta.  (1) 

To  couple  the  scattered  field  into  the  waveguide,  we  use 
Huygen’s  principle; 

(;,(r)  =  J  {pir')dGir,r')/cln 

-G{ry)df(r')ldn)dS  ,  (2) 
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where  p  is  the  nearfield  scattered  pressure  obtained  from  f 
[4].  The  calculated  pressure  field  is  represented  by  the  Green’s 
function  C.  If  we  compute  G  by  normal  modes,  we  can  write; 

G  =  (il2(I)  (Yh  ^-)  >  0) 


X,  are  the  horizontal  wave  numbers,  V|l,  are  the  modal 
amplitudes  from  which  the  pressure  field  can  be  computed, 

and  ito  =  ^  • 


where  Hq  is  an  outgoing  cylindrical  Hankel  function  and  ({)„, 
Yn  and  K„  arc  the  normal  mode  functions  and  the  vertical  and 
horizontal  eigenvalues,  respectively. 

Before  we  leave  this  section,  we  should  discuss  the  method 
we  use  for  frequency  integration.  The  inverse  Fourier  transform 
over  a  complex  field  F((o)  can  be  written; 


If  we  let  the  perturbation  -q{z)  be  added,  we  have 


d 

dz 


k  {z]-Xi  f/,  =  () 


(K) 


F(w)c‘“'r/to  . 


(4) 


We  could  use  an  FFT,  or  a  traditional  numerical  integration 
scheme,  but  we  will  present  here  a  way  that  is  more  suited  to 
our  needs  in  scattering.  For  some  small  Ato,  we  can  write 
/(r)  to  a  good  approximation  as. 


where  t/,  arc  the  perturbed  normal  modes,  and  k"(z)  = 
—q{z).  Using  completeness,  we  assume  the  following 
solution  form; 

N 

7  =  1 


/(f)  =  _L  £  I  F((o)c‘“'rf(0  .  (5) 

271  n 

If  we  can  assume  that  F((o)  doesn’t  change  much  over  any 
particular  interval,  we  can  bring  F(co,)out  of  the  integrsU 
and  analytically  evaluate  what  is  left  to  get 

71  n  ^ 


We  index  Ao)  to  allow  for  the  possibility  of  a  nonuniform 
grid  spacing.  A  similar  formulation  was  published  by 
Bracewcll  (6).  We  have  added  an  interpolation  (Newton’s 
method)  step  as  well,  although  it  is  not  always  needed.  This 
formulation  is  useful  since  we  can  accurately  compute  cases 
where  the  pressure  field  exhibits  resonance  effects  by 
concenualing  the  grid  in  their  vicinity. 


3.  NEW  PERTURBATION  APPROACH 
TO  NORMAL  MODES 

The  use  of  normal  modes  in  ocean  acoustics  is  extensively 
documented.  See,  for  example,  Boyles  [7i.  The  heart  of  the 
problem  is  to  .solve  the  depth  equation.  The  unperturbed  depth 
equation  can  be  written  as 


dz^ 


\\fi  + 


vjr,  =  0 ,  where 


(7) 


In  [2]  it  is  shown  that  this  leads  to  the  following  expression 
for  the  eigenvalues; 

X/'  —  X,’  —  </  ^  ij ,  (10) 

'■  *J 


-  an  f  * 

where  aq  s  -J-  ,  and  </,y  s  I  \)/,  q(z)  ipj  dz 

^ii  j 


(11) 


To  obtain  the  new  normal  nuxles  (/,,  we  solve  the  following 
linear  system; 

N 

ii  ^  ^  ..I 

aik  ’ 


(YX,  (X/i  +  --  ^  (/jt  <J|y  —  t/it  ,  (12) 


where  k  =  1 ,2,  ■  ■ ,  N:  i^k  , 
^ik  ~  qiJi  “  ’ 

AX,  =  Xi-X”, 

2  ,2  Jo) 

and  a,  =/t()-X,  . 
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Nowhere  do  these  expressions  involve  an  approximation  by 
tnincating  an  expansion,  as  {l(x:.s  traditional  pcrturbiition  theory. 
Accuracy  benchmark  runs  vs  SIJPFRSNAP  so  far  indicate 
near  perfect  reconsUuction  of  the  modes  for  the  cases  we 
have  tried.  Although  we  haven’t  benchmarked  computer  run 
time  yet,  our  subjective  impression  is  that  the  algorithm  is 
very  fast.  F-or  the  unperturbed  case,  we  use  an  isovehvify 
profile  over  a  half  space.  These  normal  nuKles  are  just  sine 
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functions,  and  the  q-  integrals  can  be  done  analytically. 
Bottom  attenuation  is  included  in  the  usual  way,  and  we  have 
plans  to  add  1)  sediment  layers,  2)  an  elastic  bottom,  and 
3)  extend  the  program  to  coupled  modes. 


incident  field  on  to  a  spherical  representation.  We  presented 
example  calculations  with  this  program  to  illustrate  how  we 
can  decompose  the  time  domain  response  according  to  modal 
arrivals. 


4.  DISCUSSION  OF  FIGURES 


6.  ACKNOWLEDGMENTS 


To  illustrate  the  potential  for  synthesizing  time  domain 
solutions  using  the  new  method,  we  ran  the  program  for  the 
following  case.  We  used  a  bilinear  sound  speed  profile  with 
z,  c  values  of  (0  m,  1510  m/s),  (100  m,  1476  m/s),  and  (300  m, 
1520  m/s).  The  water  and  sediment  densities  used  were  1000 
kg/m  and  1600  kg/m^,  and  the  sediment  velocity  was 
16(X)  m/s.  We  used  a  frequency  band  of  10  to  150  Hz  with 
a  A  =  0.5  Hz.  The  source  and  receiver  depths  were  both 
100  m  and  no  bottom  attenuation  was  included. 

In  Figure  1,  we  show  the  first  arrival  results  for  the  above 
waveguide.  Figure  1(a)  is  for  range  R  =  10  km,  and  (b)  is  for 
R  =  20  km.  This  shows  the  geometric  dispersion  effect  in 
waveguides.  Figure  2  (a-d)  shows  partial  modal  constructions 
for  R  =  10  km.  Figure  2(a)  is  the  time  domain  response  for 
the  single  mode  =  1,  and  2(b)  and  (2c)  are  the  single  modes 
2  and  3  respectively.  Figure  (4d)  is  the  sum  of  the  3  modes 
and  should  be  compared  to  Figure  1(a).  Figure  3(a  and  b)  are 
the  phase  and  group  velocities  calculated  for  the  first  three 
modes. 


5.  CONCLUSIONS 

We  have  presented  all  of  the  basic  features  of  performing 
pulse  scattering  in  a  waveguide  from  elastic  targets.  Primary 
emphasis  was  on  a  new  perturbation  method  for  computing 
normal  modes.  This  was  useful  in  two  ways:  1)  The  method 
is  rapid,  facilitating  the  calculation  of  many  frequencies  for 
synthesizing  a  pulse.  2)  The  perturbed  solution  can  be  written 
in  terms  of  sine  functions,  facilitating  the  projection  of  the 
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Figure  I.  (a)  Time  domain  first  arrival  solution  for  10-150  Hz,  R  =  10  km,  bilinear  profile,  (b)  Same  as  (a),  but  R  =  20  km. 
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